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Abstract 

Derivations  are  given  for  the  calculation  of  section  camberline  from  a  given  segmented  chordwise 
load  distribution,  and  for  the  inverse  problem.  In  both  cases,  the  algorithms  are  exact  within 
the  limits  of  thin-wing  theory. 

From  these  derivations,  codes  were  developed  which  are  included  in  the  section  design 
procedure  at  DREA.  The  results  from  example  cases  show  very  good  agreement  with  analytical 
solutions,  even  for  the  extreme  situation  of  infinite  leading  or  trailing  edge  slopes. 


Resume 

Presentation  de  derivations  pour  le  calcul  d’une  ligne  de  courbure  a  partir  de  segments  formes 
par  des  cordes  correspondant  a  une  repartition  de  charges,  ainsi  que  pour  le  probleme  inverse. 
Dans  les  deux  cas  les  algorithmes  sont  exacts  dans  les  limites  de  la  theorie  des  ailes  minces. 

A  partir  de  ces  derivations,  on  a  developpe  des  codes  qui  sont  inclus  dans  la  methode  de  calcul 
des  sections  du  CRDA.  Les  resultats  des  cas  donnes  en  exemple  montrent  que  la  concordance 
avec  les  solutions  analytiques  est  bonne,  meme  pour  la  situation  extreme  des  pentes  infinies  des 
bords  d’attaque  et  des  bords  de  fuite. 
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Nomenclature 


Cl  section  lift  coefficient 

cii  ideal,  or  design  lift  coefficient 

X  chordwise  coordinate/chord  (abscissa) 

y  camberline  offset,  or  normal  coordinate/chord  (ordinate) 

a  angle  of  attack 

Oj  ideal  angle  of  attack 

ao  zero-lift  angle  of  attack 

7  nondimensional  loading,  or  circulation  per  unit  span 

6  polar  coordinate 

^  chordwise  coordinate/chord 


iv 


1  Introduction 


This  document  is  concerned  with  two-dimensional,  thin,  lifting  sections  and  gives  the  derivation 
of  the  camberline  from  a  segmented  chordwise  loading  (the  design  problem)  and  the  calcula¬ 
tion  of  loading  from  a  discretely-defined  camberline  (the  analysis  problem).  To  the  author’s 
knowledge,  these  derivations  have  not  been  published  for  the  general  case,  although  the  results 
for  some  specific  ceises  are  well  known Good  examples  are  the  NACA  a-series  camberlines^, 
obtained  from  simple  two-segment  load  distributions,  which  have  been  widely  used  for  many 
yews.  /  '  / 

Modern  section  design  methods,  such  as  the  Eppler  code®,  tend  to  generate  complete  sections 
rather  than  treat  thickness  and  camber  separately.  DREA  employs  the  AADP9  code^  which  can 
be  used  in  either  fashion,  or  in  a  hybrid  of  the  two;  Reference  5  illustrates  its  use  in  a  typical 
design  exercise.  The  algorithms  presented  here  were  required  to  extend  that  previous  work 
to  the  development  of  sections  for  which  the  load  distribution  is  controlled  by  defining  the 
camberline. 

Based  on  the  present  derivations,  two  utility  programs,  CAMBER  (design)  and  CLOAD  (analysis) 
were  written  for  use  w’ith  AADP9.  Examples  using  these  programs  are  presented  and  compared 
with  analytical  results. 


2  Deriving  Camber  from  Chordwise  Loading 

This  memorandum  uses  the  standard  conventions  for  section  definition.  Thus,  the  camberline 
runs  from  the  leading  edge  at  the  origin  of  the  axis  system  to  the  trailing  edge  at  x,  or  equal 
to  1.  Camberline  offsets  are  denoted  by  y  and  lengths  are  nondimensionalized  by  the  section 
chord. 


In  thin-wing  theory®,  the  nondimensional  loading,  7(1),  on  a  thin  section  is  identified  with 
the  difference  in  local  velocity  between  the  upper  and  lower  surfaces.  From  this  definition  it 
follows  that  the  lift  coefificient  is:  ^ 

Ci  =  2  /  '){x)dx 

Jo 


If  7(1)  is  given  for  an  arbitrary  angle  of  attack,  q,  the  camberline  slope  is  found  from  the 
following  relationship®: 

In  practice,  7(2:)  is  invariably  defined  for  the  ideal  angle  of  attack,  a,,  otherwise  the  leading 
edge  loading,  7(0),  becomes  infinite.  At  the  ideal  angle  of  attack,  the  lift  coefficient  is  called 
the  section  design  lift  coefficient,  cu. 

For  this  analysis,  the  chordwise  loading,  7(1),  is  defined  in  N  segments,  so  that  the  n''* 
segment  lies  between  ordinates  i„_i  and  i„.  Therefore,  the  first  segment,  starting  at  the 
leading  edge,  lies  between  xq  =  0  and  xi,  and  the  last  segment,  ending  on  the  trailing  edge,  lies 
between  ijv-i  and  ijv  =  1,  Figure  1.  The  camberline  slope  is  then: 


7(z)di 

T 


1 


(1) 


We  proceed  by  obtaining  the  general  result  for  a  single  one  of  these  integrals  so  that; 

N 

2it 

where 


»'(f)  -  “  = E '»(«) 


n=:l 


/„(0  = /*■  ^ 


']{x)dx 

J 

The  loading  is  chosen  to  be  a  linear  function  of  x  on  each  segment,  so  that  using  the  notation 
7„  for  the  loading  at  ordinate  the  loading  on  the  segment  is: 

/  ^  iln-  7n-l)  ,  ^ 

7(a:)  =  7n  -  7- - ;; - 7  *  (^r,  -  x) 

i2:n  Xn-l) 

That  is, 

'/(l)  =  a„  +  bnX  (Xn-l  <  X  <  Xn) 

where  a„  and  are  respectively  the  loading  segment  offset  and  slope  given  by: 

an  =  ln-  and  =  7J,  (2) 

To  obtain  the  camberline  slope,  we  substitute  for  7(1)  in  the  integrals 

r  r*"  (^x^  I  /**"  xdx 

These  are  standard  integrals  and  lead  to  the  same  result  whether  ^  is  inside  or  outside  the 
segment,  that  is: 


“  ^n-l)  +  (On  +  KO  0^  [in  -  ^1  “  In  |$  -  Ir.-ij} 


(3) 


The  camberline  is  obtained  by  integrating  the  equation  for  camberline  slope.  Reversing  the 
order  of  summation  and  integration  this  gives: 

N 


or, 


where 


y(0  -  ^  S  /  In{x)d2 

-'0 

y(0  -  ^  X) 

n=l 

JniO  =  f  lrx{x)dx 
Jo 


(4) 


Therefore,  substituting  for  /n(i)i 


^(0  =  b„^{Xn  -  Xn-l) 

■f 


+a„  /  {In  |i„  -  i|  -  In  |i  - 
Jo 

+f>n  {iln|i„  -  z]  -  iln|i  -  i„_i|}di 
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Evaluating  the  integrals,  we  finally  obtain: 


MO  -  bnOr-n  -  *n-l)  (5) 

-(an  +  bnXn){^n  “  0  0^  \^n  -  ^1  -  1} 

+  (a„  +  bnXn)Xn  {In  li„l  -  1} 

-(an  +  bnXn-l){^  -  In-i)  (In  -  X„_i|  -  1} 

-(an  +  6„Xn-i)in-i  {In  -  1} 

-in— — 

For  a  =  a,,  the  section  design  lift  coefficient,  c;,,  is  obtained  directly.  Since,  for  the 
segmented  loading  (using  the  superscript  t  to  denote  that  7  is  defined  at  the  angle  of  attack), 

^  r*n 

=  2  23  / 

n=i 

simple  substitution 

cii  =  23  {2a‘„(in  -  x„_i)  +  b'„{xl  -  (6) 


In  summary,  the  camberline  slope  is  given  by  equations  (1),  (2)  and  (5),  the  cauibcrline  by 
equations  (4)  and  (5),  and  the  design  lift  by  equations  (2)  and  (6). 

It  is  easy  to  verify  that  these  results  are  correct  for  simple  cases.  For  example,  the  NACA 
(o  =  1)  mean  line  is  defined  as  corresponding  to  a  uniform  load,  7,,  at  the  ideal  angle  of  attack 
which,  since  the  loading  is  symmetrical  about  midchord,  is  zero.  Therefore, 

a\  =  7.  ,  6\  =  0  ,  zo  =  0  and  ii  =  1 

so  that 

cii  =  27. 


»’(£)  =  g{lo|l-£|- In  l£|} 

and 

v(0  =  -S{(i-0»"!i-ei  +  finie!} 
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which  is  the  result  given  by  Abbott  and  von  Doenhoff',  and  also  by  Riegels^,  who  calls  this  the 
S5  camberline. 

Finally,  in  the  general  case,  thin-wing  theory  gives  the  lift  at  zero  incidence 


Ci(a=0)  =  Cj  -  2^0 

and  the  zero  lift  angle: 

Cl 

ao  =  a~-- 

It  should  be  noted  that  if  7  is  non-zero  at  the  leading  or  trailing  edges,  then  the  camberline 
slope  at  these  points  is  infinite.  Furthermore,  the  Kutta  condition  requires  that  the  loading  go 
to  zero  at  the  trailing  edge.  In  practice,  these  strictures  are  ignored:  for  example  the  NACA 
o-series  mean  lines  are  all  loaded  at  the  leading  edge,  and  the  (a  =  1)  mean  line  is  also  loaded 
at  the  trailing  edge.  Difficulties  do  not  arise  because,  invariably,  the  leist  few  percent  of  chord 
at  each  edge  are  modified  for  manufacturing  convenience  and,  in  addition,  the  effective  camber 
at  the  trailing  edge  is  modified  by  the  boundary  layer. 


3  Deriving  Chordwise  Loading  from  Camber 


Suppose  the  loading  is  distributed  in  segments  as  before.  The  end-points  of  these  segments 
may  have  the  same  abscissae  as  the  points  defining  the  camberline.  Figure  2.  In  the  previous 
section,  we  obtained  the  following  result  for  the  camberline: 

V(0 

n=l 

Jn{i)  is  a  continuous  function  of  but  considering  only  the  discrete  values  of  ^  at  which  the 
camberline  is  defined,  we  can  replace  the  function  Jn[i)  by  discrete  values  Jm,n-  Then  the 
equation  for  the  point  on  the  camberline  is: 


1  ^ 

Vm  ~  Ch^rn  ^  ~  Jm,n 

n=l 


Since  Jm,n  can  be  expanded  in  7„  this  can  be  written  as: 


j  N-\ 

Vm  ~  ^  ^  C'f7i,n7r» 

n=0 


where  „  are  coefificients  of  the  expansion  and  the  limits  of  the  summation  have  been  shifted 
to  accomodate  the  leading  edge  (n  =  0)  and  to  exclude  the  trailing  edge  (since  the  Kutta 
condition  requires  that  7jv  =  0). 

It  is  straightforward  algebraic  manipulation  to  obteiin  Cm.n  by  evaluating  the  previously- 
given  expression  for  ^n(0.  substituting  for  On  and  and  setting  Xm  =  ^  The  final  result 
is: 


_  / 


+  E, 


m,n 


(n  =  0) 
(otherwise) 


(7) 
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where 


^n— l){ln|^Tn  ^n— ll  0 
-x„-l  {Injln-ll  -  1} 

(l„-Im)*  f,_,_  _  I  l\ 

^"2*1  -"'in'-!)  {'"I*---'-'!  -  0 
■2(x„-’x„-, ‘1 


£’m,n  =  Im  (9) 

“  (ifi  ~  ^m}  {1^  \  ^n  ~  ^m|  0 

+  I„  {In  jlnl  -  1} 

■2{r„'?^rn  f'" 

{Xm  -  Xn-l)^  /inl^  _  r  I  - 

'2(x„-x„-,)  il"l»».-»-.|-5) 

+2(x, 

Thus,  since  yo  =  0,  the  loading  7n,(o<n<Ar- 1)  is  given  by  the  solution  to  the  set  of  simulta¬ 
neous  equations; 


1 

Vm  ~~  —  2jj.  ^  -  ^m,n7n 

r»=0 


(1  <  T7l  <  A^) 


where  the  coeP  ients  Cm,n  are  defined  in  equations  (7),  (8)  and  (9). 

The  loadin,  the  leading  edge  is  finite,  according  to  thin-wing  theory,  only  at  the  ideal 
angle  of  attack,  .  the  camberline  is  invariably  defined  for  q  =  0.  If  it  is  required,  an 
approximation  for  can  be  obtained  from  the  definition®: 


where 


a.  =  -  ry[6)d9 

X  Jo 


cos  9  —  \  -  lx 
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We  replace  the  integral  by  a  summation  and  discretize  the  terms. 

a,  =  —  {arccosfl  -  2zr„)  -  arccos(l  -  2x„-i)} 

Jr  (i„  -  i„_i) 

This  approximation  may  not  be  sufficiently  accurate  for  cases  with  only  a  few  points  defining 
the  camberline  in  the  leading  or  trailing  edge  regions. 

4  Examples 

The  first  example  heis  already  been  discussed  in  Section  2:  the  constant  load  NACA  (a  —  1} 
mean  line.  The  camberline  and  slope  for  a  loading  of  7.  =  0.5  (c;,  =  1)  obtained  with  the 
CAMBER  code  were  identical  with  the  analytical  solutions  to  within  computer  precision  except 
at  the  leading  and  trailing  edges,  where  theoretically  infinite  slopes  were  encountered.  This 
situation  results  in  an  attempted  evaluation  of  the  natural  logarithm  of  zero.  In  CAMBER,  the 
problem  is  avoided  by  setting  a  lower  limit  of  10"*°  to  the  argument  of  the  natural  logarithm. 

The  next  example  is  for  a  loading  function  defined  so  that: 


n 

In 

7n 

In/lMAJC 

0 

0.00 

0.080247 

0.2 

1 

0.10 

0.320988 

0.8 

2 

0.25 

0.401235 

1.0 

3 

0.70 

0.401235 

1.0 

4 

0.90 

0.200617 

0.5 

5 

1.00 

0.000000 

0.0 

The  load  distribution  and  design  lift  were  arbitrarily  chosen  to  represent  a  typical  mean  line 
with  reduced  leading  edge  loading.  CAMBER  allows  the  input  load  distribution  to  be  normalized 
in  order  \,o  obtain  a  specific  lift  coefficient.  In  this  example,  the  values  in  the  last  column 
of  the  above  table  were  input  and  scaled  by  CAMBER  for  C(  =  0.65;  the  values  tabulated  in 
the  third  column  are  the  results.  Since  7(0)  is  finite,  these  represent  the  loading  at  the  ideal 
angle  of  attack.  The  predicted  camberline  was  subsequently  input  to  CLOAD  to  determine  how 
well  the  input  load  distribution  was  recovered.  Table  I  and  Figure  3  show  the  results  of  these 
calculations. 

There  Weis  a  small  deficit  in  loading  concentrated  at  the  leading  edge.  Other  examples  show 
that  the  effect  only  occurs  to  this  extent  when  the  loading  does  not  go  to  zero  at  the  leading 
(as  in  this  case)  or  trailing  edge.  This  is  a  numerical  artifact,  occurring  in  the  matrix  inversion 
required  to  solve  the  system  of  simultaneous  equations  (Equation  10).  It  is  beyond  the  scope 
of  this  document  to  discuss  numerical  problems  in  detail,  but  it  is  pertinent  to  note  that  a 
standard  partial  pivoting  routine  is  used  in  CLOAD,  and  that  better  methods  are  available. 

Generally,  in  this  example,  the  predicted  loading  from  CLOAD  agreed  with  the  input  to 
CAMBER  to  within  one  percent  of  the  maximum  value. 
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The  last  example  illustrates  the  evaluation  of  a  Riegels  Sj  (Birnbaum-Glauert)  camberline 
with  CLOAD.  This  is  a  cIms  of  camberlines  synthesised  from  circular  arcs  and  reflex  profiles  with 
the  general  form: 

For  this  case,  the  parameters  were  chosen  to  be:  /  =  0.1,  £i  =  -2.0  and  £2  =  0-5,  giving  a 
moderately  cambered,  highly  reflexed  camberline.  The  predictions  of  CLOAD  are  compared  with 
Riegels’  analytical  results^  in  Table  II  and  on  Figure  4. 

This  time,  since  the  loading  was  zero  at  both  leading  and  trailing  edges,  the  accuracy 
achieved  by  CLOAD  was  improved.  Agreement  between  predicted  and  analytical  loading  was 
generally  within  much  less  than  one  percent  of  the  maximum  value. 

5  Concluding  Remarks 

This  document  has  presented  derivations  of  the  relationships  between  camber  and  chordwise 
loading  for  the  general  case  in  which  these  quantities  are  expressed  at  a  finite  number  of 
discrete  points.  For  segmented  load  distributions,  typified  by  the  NACA  n-series  meanlines, 
the  derivations  are  exact  in  terms  of  thin  wing  theory.  For  most  other  cases,  they  give  very 
good  approximations,  provided  the  discretization  is  done  over  small  enough  intervals. 

Worked  examples  using  computer  codes  based  on  these  derivations  showed  very  satisfactory 
results  for  both  the  design  and  analysis  problems.  In  the  former  case,  the  numerical  imple¬ 
mentation  is  straightforward,  and  only  limited  by  computer  precision.  For  analysis,  the  present 
implementation  is  perfectly  adequate  for,  for  example,  designing  meanlines  for  propeller  blade 
sections,  but  a  robust,  high  precision  linear  solver  may  be  required  for  extremely  accurate 
results. 
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TABLE  I. 


ARBITRARY  LOADING  EXAMPLE 


Input 

CAMBER 

CLOAD 

X 

7 

y 

y' 

7 

0.0000 

0.080247 

0.000000 

0.737107 

-0.041351 

0.0050 

0.092284 

0.001168 

0.224251 

0.087432 

0.0075 

0.098303 

0.001725 

0.221690 

0.078498 

0.0125 

0.110340 

0.002827 

0.219440 

0.095316 

0.0250 

0.140432 

0.005558 

0.217869 

0.123398 

0.0500 

0.200618 

0.010973 

0.214717 

0.192775 

0.0750 

0.260803 

0.016247 

0.206002 

0.254378 

0.1000 

0.320988 

0.021184 

0.184839 

0.315498 

0.1500 

0.347737 

0.029323 

0.146886 

0.343371 

0.2000 

0.374486 

0.036035 

0.122036 

0.370805 

0.2500 

0.401235 

0.041499 

0.094693 

0.398039 

0.3000 

0.401235 

0.045526 

0.068639 

0.398409 

0.3500 

0.401235 

0.048451 

0.048879 

0.398719 

0.4000 

0.401235 

0.050448 

0.031222 

0.398968 

0.4500 

0.401235 

0.051589 

0.014502 

0.399192 

0.5000 

0.401235 

0.051901 

-0.002023 

0.399381 

0.5500 

0.401235 

0.051378 

-0.019033 

0.399568 

0.6000 

0.401235 

0.049976 

-0.037387 

0.399719 

0.6500 

0.401235 

0.047592 

-0.058661 

0.399883 

0.7000 

0.401235 

0.043975 

-0.089828 

0.400022 

0.7500 

0.351081 

0.038627 

-0.119926 

0.350013 

0.8000 

0.300926 

0.032151 

-0.138104 

0.300004 

0.8500 

0.250772 

0.024889 

-0.151944 

0.249995 

0.9000 

0.200617 

0.016942 

-0.168446 

0.200005 

0.9500 

0.100309 

0.008218 

-0.174343 

0.099892 

0.9750 

0.050154 

0.003946 

-0.166203 

0.049851 

1.0000 

0.000000 

0.000000 

-0.145228 

0.000000 

Input  7  is  interpolated  from  the  data  on  page  6;  input  c/  =  0.65 


CAMBER  output; 

Cl  =  0.650000  , 

Oi  =  0.4897  deg. 

and 

ao  =  -5.4348  deg. 

CLOAD  output: 
cii  =  0.644172  , 

Of  =  0.4393  deg. 

and 

ao  =  -5.4349  deg. 

S 


TABLE  II. 


TYPICAL  BIRNBAUM-GLAUERT  CAMBERLINE 


Input 

Analytical 

CLOAD 

z 

y 

7 

7 

0.0000 

0.000000 

0.000000 

0.028792 

0.0050 

0.001970 

0.167450 

0.173618 

0.0075 

0.002933 

0.203708 

0.202769 

0.0125 

0.004814 

0.259452 

0.260255 

0.0250 

0.009266 

0.354560 

0.358425 

0.0500 

0.017124 

0.467274 

0.470707 

0.0750 

0.023666 

0.531787 

0.533524 

0.1000 

0.028980 

0.568800 

0.572114 

0.1500 

0.036274 

0.591310 

0.595202 

0.2000 

0.039680 

0.569600 

0.572151 

0.2500 

0.039844 

0.519615 

0.521716 

0.3000 

0.037380 

0.450925 

0.452642 

0.3500 

0.032874 

0.370128 

0.371667 

0.4000 

0.026880 

0.282181 

0.283540 

0.4500 

0.019924 

0.191038 

0.192347 

0.5000 

0.012500 

0.100000 

0.101233 

0.5500 

0.005074 

0.011940 

0.013207 

0.6000 

-0.001920 

-0.070545 

-0.069288 

0.6500 

-0.008076 

-0.144999 

-0.143645 

0.7000 

-0.013020 

-0.208965 

-0.207563 

0.7500 

-0.016406 

-0.259808 

-0.258252 

0.8000 

-0.017920 

-0.294400 

-0.292766 

0.8500 

-0.017276 

-0.308510 

-0.306632 

0.9000 

-0.014220 

-0.295200 

-0.294158 

0.9500 

-0.008526 

-0.238868 

-0.234843 

0.9750 

-0.004628 

-0.179388 

-0.184229 

1.0000 

0.000000 

0.000000 

0.000000 

Analytical  results: 
cu  =  0.235619  , 

a,  =  4.2972  deg. 

and 

ao  =  2.1486  deg. 

CLOAD  output: 
cii  =  0.240397  , 

a,  =  4.2379  deg. 

and 

ao  =  2.0458  deg. 
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